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A  FORTRAN  IV  PROGRAM  TO  SCALE  AND  PLOT 
LOGARITHMIC  AND  LINEAR  DATA 


Introduction 


The  need  to  represent  data  in  a  form  more  meaningful  than  tabular  has 
prompted  the  development  of  new  software  for  general  scientific  use.  Our 
specific  needs  in  chemical  experimentation  require  the  capability  for  plot¬ 
ting  of  a  wide  variety  of  X  vs.  Y  data  types,  with  and  without  "error  bar" 
-type  error  limits  on  each  data  point.  In  addition,  it  is  often  necessary 
to  do  least  squares  calculations  on  these  same  data  and  plot  a  least  squares 
regression  line  on  the  same  plot  for  visual  indication  of  line  rity,  scatter 
and  goodness-of-fit.  Finally,  the  capability  to  provide  journal-ready  plots 
to  eliminate  the  need  for  the  user  to  take  decisions  about  scaling  is  highly 
desirable. 

Program  LPLOT  satisfies  these  basic  requirements  and  contains  such 
additional  capabilities  as  log,  semi- log  or  linear  axes,  multiple  data  sets 
on  a  single  graph,  multi-colored  plots  and  variable  origin  starting  location. 
The  main  routine  does  the  plotting;  five  subroutines  handle  the  log  and 
linear  scaling,  log  and  linear  axis  drawing  and  calculation  of  a  least 
squares  regression.  All  six  routines  are  written  in  Fortran  IV.  ^  This 
software  was  developed  to  provide  maximum  flexibility  with  a  minimum  of 
effort  by  the  user.  It  is  designed  to  run  on  a  Hewlett-Packard  1000  com¬ 
puter  system  under  an  RTE  IV  operating  system  and  makes  use  of  the  device¬ 
independent  features  of  Hewlett-Packard's  Graphics-1000  Software  as  well  as 
some  H.P.  extensions  to  standard  Fortran  IV. 

The  most  desirable  feature  of  program  LPLOT  is  that  the  calls  to  ex¬ 
ternal  routines  are  only  to  control  plotter  pen  motion,  this  makes  the  pro¬ 
gram  easily  adaptable  to  run  with  any  graphic  package. 

This  document  is  designed  to  serve  the  purposes  of  a  users'  guide  and 
an  operations  manual,  and  to  provide  sufficient  documentation  for  program 
maintenance . 

Manuscript  submitted  January  19,  1981. 
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Features  of  Program 

1.  Automatic  linear  and  logarithmic  scaling. 

2.  Choice  of  linear,  semi- log,  or  logarithmic  axes. 

3.  Least  squares  line  with  slope,  Y-intercept  and  standard  deviations 
of  each. 

4 .  Error  bars . 

5.  Multiple  plots  with  different  symbols  on  same  graph. 

6.  Data  source  can  be  from  disc  file,  cartridge  or  paper  tape  or  typed  in 
from  the  keyboard  of  a  terminal. 

7.  Output  can  be  on  graphics  terminal  or  plotter. 

8.  Plots  can  be  line,  symbols,  or  symbols  connected  by  line. 

9.  Axis  labeling  and  title  are  entered  from  terminal. 

10.  Origin  can  be  determined  automatically  or  forced  to  start  at  (0,0). 

11.  Tick  marks  are  labeled. 

12.  Multiple  color  plots  with  automatic  pen  changing. 

13.  A  manual  scaling  option  can  plot  multiple  data  sets,  which  have 
different  maximum  and  minimum  values,  on  a  single  set  of  axes. 

14.  The  routine  can  handle  up  to  256  points  per  data  set. 
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Computer  and  System  Configuration 

This  program  could  be  modified  to  run  on  many  different  computers  with 
different  plot  packages?  however,  this  routine  was  written  for  the  Hewlett- 
Packard  HP  1000  family  of  computers  using  Graphics-1000  (HP  92840A  graphics 
plotting  software) .  The  operating  system  is  RTE  IVB  with  updated  software 
revision  code  2001.  A  graphics  device,  HP  2648A  graphics  terminal  and/or  HP 
987 2B  plotter,  is  needed  to  do  the  plotting.  Older  versions  of  the  RTE  oper¬ 
ating  system  and  earlier  plotters  can  be  used  so  long  as  graphics  1000  lim¬ 
itations  are  satisfied.  EMA  and  spooling  features  were  not  utilized.  This 
routine  requires  a  26K-word  partition  to  run. 

User  Changes 

This  program  was  tested  on  a  system  probably  configured  differently 
from  that  of  the  user.  It  is  the  responsibility  of  the  user  to  make  changes 
to  logical  unit  assignments  to  implement  this  program  on  his  system.  Changes 
will  need  to  be  made  to  lines  39  through  44,  49,  299  and  390  of  program 
LPLOT.  LUG  is  the  logical  unit  number  of  the  graphics  device.  ID  is  the 
identification  number  assigned  by  the  device  link  table  to  a  graphics 
logical  device.  Line  68  of  program  LPLOT  does  not  permit  a  logical  unit 
number  greater  than  30.  This  limit  was  set  to  avoid  input  errors  and  may 
need  to  be  changed  by  the  user. 

Program  Background 

The  Hewlett-Packard  Graphics  1000  plot  package  (HP  92840A)  does  not 
provide  logarithmic  plotting,  linear  or  logarithmic  scaling.  It  is  left 
for  the  user  to  determine  the  limits  of  the  data  and  how  the  axes  of  a  plot 
are  to  be  annotated. 

The  existing  graphics  software  forces  the  user  to  draw  two  linear  axes 
with  one  call  to  an  axes  drawing  routine.  This  action  restricts  both 
axes  anotation  to  one  format.  Program  LPLOT  is  not  limited  by  this  restriction 
because  it  uses  two  calls,  one  for  each  axis,  so  each  axis  can  be  anotated 
differently . 
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The  need  to  represent  data  logarithmically  as  opposed  to  linearly 
is  evident  when  data  spans  many  orders  of  magnitude  or  conforms  to  logarithmic 
interpretation . 

It  is  because  of  these  needs  and  restrictions  program  LPLOT  was 
written. 

Cautions 


Erroneous  control  codes  in  a  data  file  or  failure  to  separate  data  with 
a  comma  or  space  can  produce  unpredictable  plotting  results. 

Data  files  set  up  for  error  bar  plotting  must  set  DELX  and  DELY  equal 
to  zero  for  those  data  points  where  error  bars  are  not  desired. 

The  subroutine  SCALE  transforms  the  array  A  into  exponential  notation, 
storing  the  exponent  in  IEXP.  This  process  changes  the  array  A  by  some 
factor  of  10.  If  the  user  calls  this  subroutine  from  any  other  program  he 
should  be  aware  that  any  value  passed  back  to  the  calling  routine  could  be 
changed.  It  is  for  this  reason  A  and  IEXP  must  be  used  together. 

A  symbol  type  plot  cannot  center  a  symbol  directly  on  an  axis  because  the 
plotting  window  boundary  is  set  on  the  axis  and  calls  cannot  be  made  to  plot 
outside  this  area. 

Some  errors  in  input  data  cam  generate  error  conditions  in  Hewlett- 
Packard  library  routines.  As  written,  the  program  allows  these  to  be 
printed  on  the  standard  list  device  (logical  unit  6) .  To  avoid  this,  the 
user  must  supply  his  own  error  routine  as  described  in  relocatable  library 
routine  ERQ.E.^ 

The  user  should  be  aware  that  significance  in  single  precision  real 
format  is  to  six  or  3even  decimal  digits.  It  is  due  to  this  limitation 
that  plotting  errors  will  be  encountered  if  input  data  exceedes  this  limit. 
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Input 


Data  can  be  input  to  the  program  from  a  disc  file,  cartridge  or  paper 
tape  or  typed  in  from  the  keyboard  of  a  terminal.  Disc  files  must  be  type 

3  or  4.^  The  program  can  handle  an  array  up  to  256  data  points.  The 

(4) 

format  is  X,  Y,  DELX,  DELY  using  free  field  input.  The  X  and  Y  pair  is 

the  position  of  the  point  along  the  X  and  Y  axis  respectively.  The  optional 

pair  DELX  and  DELY  is  the  standard  deviation  in  X  and  Y.  These  values  are 

used  for  drawing  error  bars  and  will  be  doubled  and  scaled  to  plotter  units 

to  provide  the  horizontal  or  vertical  separation  of  error  bars.  All  numeric 

•38  38 

input  data  must  be  in  the  range  of  10  to  10 

Loading  Procedure  on  HP  1000 

The  loader  must  be  loaded  as  a  large  background  program  (type  4)  and 
will  require  a  size  of  26  pages  to  accommodate  the  Graphi 'S-1000  routines. 
After  compiling  &LPLOT ,  &SCALE,  &LSCAL ,  &LINAX,  &LOGAX  and  &LSREG ,  execute 
the  loader  on  the  command  file  or  interactively  as  follows : 


RU.LOADR 
SZ.26 
OP,  LB 
RE, %LPLOT 
RE , % SCALE 
RE, %LSCAL 


RE, %LINAX 
RE , % LOG AX 
RE,%LSREG 
RE , %DLTBL 
SE, %GPS 
END 


where  %DLTBL  is  the  device  link  table  and  %GPS  is  the  Graphics-1000  library 
file  created  when  Graphics-1000  was  loaded. ^  Alternatively,  the  source 
version  available  from  the  author  contains  a  loader  command  file  which  can 
be  used.  The  loaded  program  will  occupy  approximately  26K  words  of  memory. 
It  would  be  possible  to  decrease  the  main  program  space  somewhat  by  using 
EMA^6'  for  array  space.  However,  the  bulk  of  this  26K  words  is  required 
for  Graphics-1000  routines. 


Program  Source  Availability 

A  program  source  is  available  from  the  author  on  a  user-supplied 
Hewlett-Packard  264X-type  cartridge  tape.  This  tape  contains  8  files.  The 
first  file  is  a  description  of  what  is  on  the  tape.  The  second  file  is  a 
command  file  which  may  be  used  by  the  loader  to  load  the  programs.  The 
third  file  is  the  source  for  program  LPLOT.  The  fourth  file  is  the  source 
for  subroutine  SCALE.  The  fifth  file  is  the  source  for  subroutine  LSCAL. 

The  sixth  file  is  the  source  for  subroutine  LINAX.  The  seventh  file  is 

the  source  for  subroutine  LOGAX.  The  eight  file  is  the  source  for  subroutine 

LSREG. 


Error  Messages 

Program  LPLOT  checks  for  errors  that  could  occur  when  reading  a  disc 
resident  file.  Messages  to  the  user  are  sent  with  the  name  of  the  file 
if  a  problem  is  encountered.  The  program  has  several  built  in  checks  to 
catch  typing  errors  by  the  user. 

Testing  of  Program 

The  author  has  used  a  wide  variety  of  data  types,  multiple  files,  and 
several  input  sources  to  debug  this  program.  In  addition,  the  program  has 
been  used  extensively  by  other  individuals  with  different  applications. 

Execution  of  Program 

The  execution  of  LPLOT  can  be  performed  by  supplying  terminal  and 
graphics  device  Logical  Unit  Numbers  (LU  and  LUG)  via  the  run  string 
parameters.  These  parameters  (globals)  are  retrieved  by  a  call  to  the 
Hewlett-Packard  routine  RMPAR  in  LPLOT.  If  these  parameters  are  not  supplied, 
the  program  retrieves  LU  and  prompts  the  user  for  the  graphics  device  LUG. 
Device  selection  is  to  be  made  between  the  plotter,  for  a  hard  copy,  or  the 
graphics  terminal.  The  data  source  is  from  a  disc  file,  cartridge  or  paper 


tape  or  entered  via  the  keyboard  of  a  terminal.  If  the  source  is  a  disc 
file  the  name  of  the  file  is  requested.  If  the  source  is  other  than  the 
disc  a  logical  unit  number  is  requested.  One  of  three  types  of  plots  can 
be  selected;  a  straight  line  connecting  each  point,  a  symbol  at  each  point 
or  a  symbol  at  each  point  with  connecting  lines.  The  user  may  select  either 
axis  to  be  linear  or  logarithmic.  Labels  are  entered  from  the  users  terminal, 
any  ASCII  character  is  permitted  (capital  and  lower  case  letters,  numbers 
and  symbols) .  The  X  and  Y  axis  labels  cannot  exceed  30  characters  and  the 
title  of  the  plot  cannot  exceed  40  characters  in  length;  characters  beyond 
these  limits  are  ignored.  The  starting  position  of  the  origin  can  be  forced 
to  start  at  X  =  0  and  Y  =  0  or  the  user  can  let  the  scale  routine  determine 
an  origin  that  will  maximize  the  size  of  the  plot  vs.  the  size  of  the  axes. 
Scaling  is  performed  automatically,  however  a  manual  override  is  provided. 

The  manual  scaling  mode  is  used  to  increase  the  limits  between  the  maximum 
and  minimum  values  of  an  axis.  This  feature  is  necessary  when  plotting 
multiple  plots  on  one  graph  when  the  maximum  and  minimum  values  of  all 
of  the  data  sets  are  not  within  one  data  set.  When  plotting  multiple  data 
sets  the  first  set  plotted  must  have  the  smallest  and  largest  values  of 
all  of  the  data  sets.  If  this  condition  cannot  be  met  the  user  must  specify 
manual  scaling  and  enter  minimum  and  maximum  values  of  the  entire  set  of 
data.  A  least  squares  line  can  be  drawn  on  the  plot  with  slope,  y-intercept 
and  respective  standard  deviations  printed  on  the  users  terminal  or  printer. 
Error  bars  can  be  drawn  around  each  point  provided  requirements  in  the  in¬ 
put  section  of  this  paper  have  been  met.  Error  bars  that  are  small  enough 
to  distort  the  symbol  printed  at  a  data  point  are  suppressed  and  a  message 
is  printed  on  the  users  terminal  of  this  action.  Six  or  less  plots  can  be 
made  using  the  same  set  of  axes  provided  all  data  sets  fall  within  the  limits 
of  the  first  data  set  plotted.  The  user  has  the  option  to  make  pen  color 
changes  when  doing  multiple  plots  on  the  same  axes. 

Subroutines 


The  axis  drawing  and  scaling  subroutines  of  program  LPLOT  do  not  make 
calls  to  external  routines  other  than  for  plotter  pen  movements. 
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The  subroutine  SCALE  uses  a  table  look-up  method,  based  on  the  difference 
between  maxinum  and  minimum  values ,  to  determine  axis  scaling  and  number  of 
tick  marks  to  be  placed  on  each  axis.  The  SCALE  routine  uses  an  algorithm 
that  shifts  decimal  points  to  increase  numbers  that  are  less  than  one  and 
decrease  numbers  that  are  greater  than  1000. 

The  subroutine  LSCAL  logarithmically  scales  an  array  to  determine  the 
exponent  to  which  the  base  10  must  be  raised  to  produce  a  number  less  than 
or  equal  to  the  minimum  value  of  the  array  and  greater  than  or  equal  the 
maximum  value  of  the  array. 

The  subroutine  LINAX  draws  a  linear  axis  and  labels  tick  marks  in  the 
X  or  Y  direction. 

The  subroutine  LOGAX  draws  a  lograthmic  axis  and  labels  tick  marks  in 
the  X  or  Y  direction. 

The  subroutine  LSREG  does  a  ieast  squares  regression  calculation  including 
the  standard  deviation  of  the  slope  and  y  intercept. 

Plotter  Set  Up 

The  HP  987 2B  plotter  is  used  when  hardcopy  results  of  LPLOT  are  re¬ 
quired.  The  plot  generated  by  LPLOT  is  designed  to  fit  on  standard  8*j  X  11 
inch  paper  with  adequate  margins  for  publication.  In  order  to  center  the 
plot  the  paper  must  be  placed  horizontally  *»  inch  from  the  left  side  and  1 
inch  from  the  bottom  of  the  plotter  bed.  This  displacement  is  necessary 
because  the  HP  9872B  plotter  with  advance  option  OFF  places  the  lower  left 
corner  of  the  plotting  window  at  (520,  380)  instead  of  (0,0).  For  multi¬ 
colored  plots  pen  placement  is  as  follows;  pen  1  is  black,  2  is  red,  3  is 
green,  and  4  is  blue.  If  the  plot  is  to  be  one  color  that  color  pen  must 
be  in  pen  holder  1. 
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Plotting  Accuracy 

The  HP  9872B  plotter  is  divided  into  plotter  units  where  one  unit  = 
0.025  mm.  This  is  to  say  that  the  overall  resolution  of  the  plot  is  one 
part  in  N  where  N  is  the  number  of  plotter  units  occupied  in  the  X  or  V  dir¬ 
ection.  The  X  axis  is  7.7  inches  (7,823  plotter  units)  and  the  Y  axis  is 
5.25  inches  (5,334  plotter  units).  The  plotter  resolution  is  one  part 
in  7,823  in  the  X  direction  and  one  part  in  5,331  in  the  Y  direction. 
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Mnemonic  List 


BARX  -  Minimum  value  of  X  error  bar 

BARY  -  Minimum  value  of  Y  error  bar 

DELX  -  Experimental  error  in  X  (standard  deviation) 

DELY  -  Experimental  error  in  Y  (standard  deviation) 

DIF  -  Difference  between  maximum  and  minimum 
HH  -  Half  height  in  character  cells 
HW  -  Half  width  in  character  cells 
IAX  -  Label  of  X  axis  (30  characters  max) 

I AXES  -  Type  of  plot  (log  and/or  linear) 

IAY  -  Label  of  Y  axis  (30  characters  max) 

XBAR  -  Type  of  error  bars 

ID  -  Identification  number 

IDONE  -  Check  for  termination  or  multiple  plot 

IEXP  -  Exponent  of  base  10  in  X  data  scale 

IFMT  -  Source  of  data 

IHED  -  Title  of  plot  (40  characters  max) 

I LINE  -  Check  for  least  squares  line 
IPEN  -  Pen  number  of  plotter 

ISCAL  -  Set  to  zero  for  automatic  scaling,  1  for  manual 

I ZERO  -  Set  to  zero  to  force  origin  to  start  at  zero 

JCHAR  -  Character  to  be  plotted  at  each  data  point 

JEXP  -  Exponent  of  base  10  in  Y  data  scale 

JJ  -  Type  of  plot 

KK  -  Number  of  plots  on  same  axes 

LU  -  Logical  unit 

LUG  -  Logical  unit  of  graphics  device 

LUT  -  Logical  unit  of  cartridge  tape  or  keyboard 

NAME  -  Name  of  data  file 

NOBAR  -  Check  to  see  if  error  bar  was  too  small  to  plot 
NP  -  Number  of  points 

51  -  Standard  deviation  of  slope 

52  -  Standard  deviation  of  Y  intercept 
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SLOPE  -  Slope  of  least  squares  line 
SXTIC  -  Interval  between  X  tick  marks 
SYTIC  -  Interval  between  Y  tick  marks 


X  -  Displacement  along  X  axis 

XBAR1  -  Distance  to  left  of  character  of  X  error  bar 
XBAR2  -  Distance  to  right  of  character  of  X  error  bar 
XEND  -  X  value  at  end  of  least  squares  line 
XMAX  -  Maximum  value  of  X 
XMIN  -  Minimum  value  of  X 

XST  -  X  value  at  start  of  least  squares  line 
XTIC  -  Number  of  tick  marks  on  X  axis 
Y  -  Displacement  along  Y  axis 
YBAR1  -  Distance  above  character  of  Y  error  bar 
YBAR2  -  Distance  below  character  of  Y  error  bar 
YEND  -  Y  value  at  end  of  least  squares  line 
YINT  -  Y  intercept  of  least  squares  line 
YMAX  -  Maximum  value  of  Y 
YMIN  -  Minimum  value  of  Y 

YST  -  Y  value  at  st3urt  of  least  squares  line 
YTIC  -  Number  of  tick  marks  on  Y  axis 
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Formulas 


The  following  formulas  were  used  in  subroutine  LSREG: 

Slope  =*  (DIF)(IXIY)  -  (IXI Y) 

(DIF) (EX) 2  -  (IX)2 


Y- Intercept 


Standard  Deviation 
of  Slope 


Standard  Deviation 
of  Y  Intercept 


DIF  ■  the  interval  over  which  the  calculation  is  computed  (FROM-TO+1) 
IX  *  summation  of  X  values 
IY  »  summation  of  Y  values 


INITIALIZE 
PLOT  PACKAGE 


start 


PROGRAM 

SCALE 
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PLOTTING  EXAMPLES 


LINEAR  AXES 

LINE  PLOT 

DATA: 

391 .51282 

.25431687 

397.39148 

.33236635 

403-27008 

. 40419948 

409. 14874 

.47188836 

415.02740 

.50006902 

420.90607 

.50377190 

426.78473 

.49136621 

432.66333 

.45517337 

438.54199 

.40295619 

444.42065 

.34763810 

450.29932 

.29423952 

456.17798 

.24892941 

462.05664 

.21922916 

467.93524 

.19685039 

473.81390 

.18386519 

479.69257 

.17585301 

485.57123 

.17474783 

491 .44989 

. 17474788 

497.32849 

. 17751071 

503.20715 

.18414146 

509.03582 

.  19436386 

514.96448 

.20721093 

520.84314 

.23456278 

526.72180 

.26453930 

532.60046 

.30777734 

538.47913 

.  35640281 

544.35767 

.41207349 

550.23633 

.47575635 

556.11499 

.53598559 

561.99365 

.59787261 

567.87231 

.65451026 

573-75098 

.69691944 

579.62964 

.72413313 

585.50330 

.74057186 

591.38696 

.73394120 

597.26550 

.71529221 

603.14417 

.67302108 

609.02283 

.61983705 

614.90149 

.56016028 

620.78015 

.49440533 

626.65881 

.42243403 

632.53748 

.35115349 

638.41614 

.28526038 

644.29480 

.23097113 

650.17346 

. 18220749 

656.05212 

. 1 Uo98l63 

661 .93079 

.12170190 

667.80933 

.978035632-0 1 

673.68799 

.63851 367E-0 1 
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o7  9 . 5oob5 
oo5 . 445 j 1 

3*1.32397 

0*7.20264 
703 . Ool 30 
703.9599o 
7  14. d  3^02 
720.71710 

2.  LINEAR  AXES 

SYMBOLS  WITH  CONNECTING  LINES  PLOT 

FOUR  DIFFERENT  SYMBOL  PLOTS  ON  SAME  AXIS 

DATA: 

1 .  2,10.1 

4.10.2 
o,10.2 
0,10.1 

10.10.2 

12,10.2 

14.10.3 

16. 10.3 

16.10.3 

20.10.3 

22.10.4 

24.10.4 
2b , 10.5 

26. 10.4 

30,10.2 
32,10.0 

34.9.94 

3b, 10.0 

38.9.95 

40.10.1 

2.  2,10.0 

4.10.1 

6.9.94 

6,10.0 

10.10.1 
12,10.1 

14.10.3 
1b, 10.2 
18,10.0 
20,  10.3 

22.10.3 

24.10.4 

2b,  10. 3 

28,10.3 

30,10.1 

32.9.91 
34,9.82 

36.9.92 
38,9.68 

40.9.94 


.716Q4985E-01 
.O5202370E-0  1 
.O0505591E-Q1 
.54427408E-01 
.519406732-01 
. 4676  3640E-0  1 
.  487 6  3640E-0  1 
.47244094E-01 


3.  2,10.3 
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4.10.6 
c  ,  10.6 
0,10.6 

10.10.9 

12.10.9 
14,11.0 
16,11.1 

13.11.1 

20.11.2 

22.11.4 
24,  11 .u 
26,11.0 

20. 1 1.5 

30.11.5 

32.11.5 

34,1  1.5 

36.11.5 
33,  1  1.5 

40.11.5 


4.  2,10.2 

4.10.2 

0,10.3 

3,10.5 
10,  10.4 

12.10.4 

14.10.4 

16.10.5 
13,10.4 
20,  10,5 
22, 10.7 

24.10.7 

26.10.8 
23,  10.6 

30.10.3 

32.10.3 

34.10.2 

36.10.2 

33,10.1 
*0, 10.3 

3.  SEMI-LOG  AXES  (LOGARITHMIC  X  .  LINEAR  Y ) 
SYMBOLS  WITH  CONNECTING  LINES  PLOT 
A  SINGLE  PLOT  DIVIDED  INTO  3  OVERLAPPING  DATA 
rloES  WiiH  A  LcA^T  SCL'ARtS  wINE  0.^  z. AC H  FIL: 
DA  i  n  • 

.  0  0  1  o  .  2  4  5  c 

.0040,2270 
.0095,2033 
.0155. 1790 

•034Q, 1=90 

.  1  175,  1535 
1 .24, 147o 
3-33,  1  uU5 
32.4, 1245 

105.3.1033 


445 . 2 , 867 
890.4, o35 
1480,465 


4.  LOGARITHMIC  AXES 

SYMBOL  PLGT  WITH  X  AND  Y  ERROR  BARS 
LEAST  SQUARES  LINE 
DATA : 

.023,124, .009,43.3 
.79,435,-26,153.5 
3.5,1255,1 .54,330 
43.8,10235,7.6,3630 
125, 120450,22.4,32540 
340, 1350450, 165,430000 
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PLOT  2  «  O-B  .  +-P  .  *-Cd 


C\I 


* 

♦ 


* 

* 


♦ 

i 


X  G  r 

‘X  Or 
X  Oi 
X  0  i 
X  0f 

.x.o>' 
x-‘  oV 
x'  o  > 

X  ($ 

x  dj 1 1 

X  $.  ’ 

X  0  r 
X  0  f  ’ 

X© 

♦  x  ©> 

*.  x  Ot 

>  x  o  't 
*  xo  > 

Efl  f’ 

'  ’*  x  6)  '+ 


a 


Q 


Q 

m 


Q 

OJ 


G 


oi 


a 


<Wdd>  NCIiVdlN33NC: 
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TIME  IN  MINUTES 


AlIAIlddOSSV  avion 


PLOT  4  <  EXAMPLE  OF  A  LOG  LOG  PLOT 


LOG  X  AXIS 


References 


1.  Hewlett-Packard,  RTE  Fortran  IV  Reference  Manual,  January  1980. 

2.  Hewlett-Packard,  RTE  Relocatable  Library  Reference  Manual,  December  1978, 
p.  3-7. 

3.  Hewlett-Packard,  HP  92068A  RTE-IVB  Terminal  User’s  Reference  Manual, 
January  1980,  p.  3-7. 

4.  Hewlett-Packard,  RTE  Fortran  IV  Reference  Manual,  January  1980,  p.  7-4. 

5.  Hewlett-Packard,  HP  92840A  Graphics  Plotting  Software  User's  Manual, 
January  1980,  p.  5-1. 

6.  Hewlett-Packard,  HP  92068A  RTE-IVB  Terminal  User's  Reference  Manual, 
January  1980,  p.  1-2,  4-56. 


SAMPLE  DIALOGUE  FOR  EXECUTION  OF  PROGRAM  LPLOT 


: RU, LPLOT 

OUTPUT  ON  GRAPHICS  TERMINAL 
ON  PLOTTER 

0 

***  SOURCE  OF  DATA  *** 
DISC  FILE  TYPE  0 

CARTRIDGE  TAPE  TYPE  1 

TERMINAL  TYPE  2 


TYPE 

TYPE 


ENTER  NAME  OF  DATA  FILE 
P  VP  1 

***  TYPE  OF  PLOT  *** 

LINE  PLOT 

SYMBOLS  CONNECTED  WITH  LINES 
SYMBOLS  PLOT 

2 

***  TYPE  OF  AXES  *** 
LINEAR  X  AND  Y 
LOGARITHMIC  X  AND  Y 
LOGARITHMIC  X  AMD  LINEAR  Y 
LOGARITHMIC  Y  AND  LINEAR  X 
3 

"enter  X-LAEZL,Y-LASEL  4  TITLE  0 

TIME-MINUTES 

LOG-ABSORBANCE 

?VP  WITH  FE+3 

***  ORIGIN  LOCATION  *** 

TO  FORCE  LINEAR  AXIS  TO  START  A 
TO  LET  SCALE  DETERMINE  ORIGIN 

1 

***  SCALING  *** 

AUTOMATIC  TYPE  0 

MANUAL  TYPE  1 

0 

LEAST  SQUARES  LINE  ? 

NO  TYPE  0 
YES  TYPE  1 


TYPE 

TYPE 

TYPE 


TYPE  0 
TYPE  1 
TYPE  2 
TYPE  3 


SEPARATE  LIN 


TYPE  0 
TYPE  1 


ERROR  BARS 


NONE 

X 

Y 

X  4  1 


TYPE 

TYPE 

TYPE 

TYPE 


TO  EXIT 

IF  YOU  WANT  ANOTHER  PLOT  ON  SAME  AXES 
AND  ALL  X  4  Y  VALUES  ARE  WITHIN  THE 
SCALES  OF  THE  FIRST  PLOT 


TYPE 


TYPE 


PAGE 

0001 

FTN.  9:10  AM  MON.,  26  JAM.,  1931 

-  ^ 

0001 

FTN4 , 

L 

0002 

PROGRAM.  LPL0T(4,99) 

0003 

C 

LOGARITHMIC  4  LINEAR  PLOTTING  ROUTINE  FOR  THE  HP  98723  PLOTTER 

0004 

C 

MAXIMUM  256  PTS.  ALLOWED,  X  AND  Y  MAY  3E  INPUT  FROM  A  DISC 

0005 

c 

FILE.  CARTRIDGE  TAPE  OR  KEYBOARD  FORMATTED  IN  (X, Y)  PAIRS 

0006 

c 

FOR  EXPERIMENTAL  ERROR  BARS  FORMAT  IS  ( X , Y , DELTAX , DELTAY ) . 

0007 

c 

0008 

COMMON  IGCB0921 

0009 

DIMENSION  IDCE ( 1 44 ) , NAME ( 3 ) , I AX ( 1 5 ) , I AY ( 1 5 ) , IHED( 20 ) , IBUF (40) 

0010 

DIMENSION  X(25o) ,Y(25o) ,DELX(256) ,DELY(256) 

0011 

DIMENSION  IP  RAM( 5 ) , IOBUF ( 20 ) , JCHAR ( 6 ) 

0012 

DATA  JCHAR/ 1  HO,  1 H+ ,  1  H!« ,  1 HX ,  1 H# ,  1H$/ 

0013 

CALL  RMPAR( IP RAM) 

0014 

LU=IPRAM( 1 ) 

0015 

LUG=IPRAM(2) 

0016 

c 

0017 

C  SET 

DEFAULT  VALUES 

0018 

C 

0019 

ILUG=0 

0020 

iF>rr=o 

0021 

JJ=1 

0022 

IAXES=0 

0023 

IZERO= 1 

0024 

ISCAL=0 

0025 

ILINE=C 

0026 

IBAR=0 

0027 

IDONE=0 

0023 

IP  e/a  7 

0029 

LERR=0 

0030 

r 

0031 

C  ESTABLISH  OUTPUT  DEVICE 

0032 

r* 

O 

0033 

IF (LUG. ME. 0) GOTO  15 

0034 

WRITE(LU, 10) 

0035 

10 

FORMAT ( "  OUTPUT  ON  GRAPHICS  TERMINAL  TYPE  0",/, 

0036 

1  "  ON  PLOTTER  TYPE  1") 

0037 

READ(LU,*)IL’JG 

0038 

C 

DEFINE  LU  AND  ID  NUMBERS  OF  PLOTTER  AND  GRAPHICS  TERMINAL 

0039 

LUG=24 

0040 

IF ( ILUG . EC . 1 ) LUG=20 

0041 

15 

KK=Q 

0042 

IF(LU.LE.O)LU= 1 

0043 

ID=1 

, 

0044 

IF ( LUG . EQ . 20 ) ID=2 

r 

0045 

C 

TO  CENTER  CHARACTER  SET  HALF  WIDTH  4  HALF  HEIGHT 

0046 

C 

DEPENDING  ON  IF  USING  PLOTTER  OR  CRT 

0047 

;iW=0.5 

0048 

HH=0 . 5 

0049 

IF(ID.EQ.I)  GOTO  20 

0050 

HW=0.333 

0051 

hK=0 . 25 

0052 

20 

KK=KKf1 

0053 

C 

i 

0054 

C  DETERMINE  DATA  SOURCE  AND  READ  DATA 

i 

I 

0055 

1 

26 


PAGE 
0056 
0057 
005d 
0059 
OOoO 
0061 
0062 
00b  3 
006*1 
0065 
OObb 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0078 
0079 
3080 
0081 
0082 
00o3 
0084 
0085 

3086 

3087 
0088 
0089 
0090 
0091 
0092 
0093 
009** 
0095 
0096 
0097 
0098 
3099 
0100 
0101 
0102 
0103 
0104 
0105 
Cl  06 
0107 
0108 
0109 
0110 


0002  LPLOT  9:10  AK  MON.,  26  JAN.,  1981 
25  'WRITE(LU,30) 

30  FORMAT ( "  ***  SOURCE  OF  DATA  ***"/ 

1  ”  DISC  FILE  TYPE  0"/ 

2  "  CARTRIDGE  TAPE  TYPb  1"/ 

3  "  TERMINAL  TYPE  2") 

READ  (LU,»)  IFWT 
IF(IFMT.EQ.O)  GOTO  70 

35  IFCIFMI.EQ. 1 )WRITS(LU,40) 

IF ( IFMT . NE . 1 ) WRITE ( LU , 45 ) 

40  FORMAT (”  ENTER  LU  NUMBER  OF  CTU”) 

45  FORMAT ( "  ENTER  LU  NUMBER  OF  TERMINAL") 

READ(LU, *)LUT 
IF(LUT.GT.30)  GOTO  35 
50  WRITE (LU, 55) 

55  FORMAT ( "  ENTER  NUJ-SER  OF  DATA  POINTS  (MAX=256)") 

READ(LU,*)MP 
IF(NP.GT.256)  GOTO  50 
IF  ( IFMT .  NE .  1 ) 'WRITE  ( LU ,  60 )  NP 
60  FORMAT ("  ENTER  ",I3,”  DATA  POINTS") 

D065  1=1, NP 

65  READ  (LOT.*)  X( I) ,Y( I) ,DELX(I) ,DELY( I) 

GOTO  110 

73  'WRITE  (LU,  75) 

75  FORMAT ("  ENTER  NAME  OF  DATA  FILE”) 

READ(LU,80)NAME 
30  FORMAT ( 3 A2) 

CALL  QPEN( IDCB , IERR , NAME , 3) 
IF(IERR.LT.0)WRITE(LU,85)NAME 
IF ( IERR. LT.O) GOTO  25 

85  FORMAT ("  FILE  ",3A2,"  DOESN'T  EXIST  OR  IS  ALREADY  OPEN" 

K=  1 

90  D095  J=1 ,40 

95  IBUF(J)=0 

CALL  READF ( IDCB , IERR , I5UF , 40 , LEN ) 

IF(LEN.EQ.-I)  GOTO  100 
IF( IERR . LT . 0 ) WRITE( LU , 1 05 ) IERR 
CALL  CODE 

READ( I3UF , * ) X ( K ) , Y ( X ) , CELX ( K ) , DELY( X ) 

X=K+1 
CO  TO  90 
100  NPrK-1 

CALL  CLOSE( IDCB, IERR) 

IF(IERF.LT.0)WRITE(LU,105)IERR 
105  FORMAT ( "  DISC  I/O  ERROR  =  "14) 

C 

C  SCALE  DATA  TO  FIRST  PLOT  (FOR  MULTIPLE  PLOTTING  ONLY) 

C 

110  IF(KK.EQ.I)  GOTO  165 
IF(  IAXES.EQ.  1  )GCTO  1**5 
IF( IEXP.EQ.O.AND.  JEXP.EQ.OGCTO  0*5 
D0 140  :<=i,np 
IF(IEXP) 115, 125, 120 
115  X(K)=X(K)/( 10.**IA£S( IEXP) ! 

GOTO  125 

’20  X(K)=X(K’,*O0.##IAB3CEXP)) 
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PAGE 

30C3 

-PLOT  9:10  A.-!  MON..  26  JAN.,  1: 

?3i 

01" 

125 

IF(JE:<P)130, 140, 125 

0112 

1<0 

Y( K)=YC<)/ ( 1C . **IA8S( JEXP ) ) 

0H3 

GOTO  1-tC 

01 14 

125 

Y(K)=Y(K)*( 10.**IABS( JEXP ) ) 

01  T 

140 

CONTINUE 

011c 

145 

IF'.IAXES.NE.  1  .AND.IAXES.NE.2)GOTO 

155 

0117 

XI 50  K=1,NP 

ana 

X( K ) = ALOGT ( X( K } ) 

0119 

IF(DELX(K) .EQ.O. )GOTO  150 

0 1 20 

DSLX ( K )  =  ALOGT ( .AES ( CELX (K) ) ) 

0121 

150 

CONTINUE 

0122 

155 

IF(IAXES.NE.1.AND.IAXES.NE.3)GOTC 

165 

0123 

1X0160  K=1,NP 

0124 

Y(K)=ALOGT(Y(K)) 

0125 

IF(DELY(K) .EQ.O. )GOTO  160 

012b 

DEL Y  ( K )  =  ALOGT  ( ABS  ( DEL Y  ( K ) ) ) 

0127 

IbO 

CONTINUE 

0123 

c 

0129 

C  DETERMINE  PLOTTING  FEATURES 

0130 

C 

0131 

165 

WRITE (LU, 170) 

I 

0132 

170 

FORMAT ( "  ***  TYPE  OF  PLOT  *«",/ 

0132 

1  "  LINE  PLOT 

TYPE  1 " , / 

0134 

2  "  SYMBOLS  CONNECTED  WITH  LINES 

TYPE  2",/ 

0135 

3  "  SYMBOLS  PLOT 

TYPE  3”,)  j 

0126 

READ^LU,*) JJ 

0127 

if(kk.gt.i;goto  300 

0126 

WRITE! LU, 175) 

0139 

175 

FORMAT ( "  ***  TYPE  OF  AXES  *»•",/ 

0140 

1  "  LINEAR  X  AND  Y 

TYPE  0",/ 

3141 

2  "  LOGARITHMIC  X  AND  Y 

TYPE  1",/ 

0142 

2  "  LOGARITHMIC  X  AND  -INEAR  Y 

TYPE  2"./ 

0142 

4  "  LOGARITHMIC  Y  «ND  LINEAR  X 

TYPE  3”) 

0144 

=EAD(LU.*  TAXES 

0145 

WRITE  (L'J,  ISO) 

014b 

180 

FORMAT ( "  ENTER  X -LABEL,  Y-LABEL  !t 

IITLE  ON  3  SEPAFATE  LINES") 

0147 

READ(LU,1S5>  IAX,I*Y,IKED 

014C 

*85 

FORMAK 15A2/.15A2/.20A2) 

0 1  uv 

IF' I AXES. EC. 1 )GOTO  ’95 

3150 

WRITE( LU, ibO) 

0151 

•4Q 

FORMAK"  •**  ORIGIN  LOCATION 

/ 

»  • 

0152 

1  ”  TO  FORCE  -INEAP  AXIS  TO  START  AT  0  TYPE  0",/ 

3153 

2  "  TO  LET  SCALE  DETERMINE  ORIGIN 

TYPE  1") 

015- 

READ( LU , * ) IZERC 

0155 

195 

WRITE! LC , 200 I 

0156 

2CG 

FORMAK”  *»•  SCALING  *«•",/ 

oi57  i  "  automat ic  type  O",/ 

0150  2  "  ’1ANUAL  TYPE  1") 

01;:*  PEADfLU,*'  IC'AL 

0 1  ~0  IF '  I5CAL .  EC  .  )'•  GOTO  ,  ’  0 

0T1  WRITE! -L,. >05 

0‘-2  205  FCRNA!"’  tNTER  :KIN,>KAX,YM.:N,YNAX") 

I '  4  3  REAL < L'J ,  • ' XMIN .  /MAX .  YMN ,  YVAX 

01-.4  : 

0155  0  INITIALIZE  PLOTTER 
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PAGE  0004  LPLG1  y:lD  W  KGN . ,  26  JAI..,  1981 
0166  0 

01 67  £10  CALL  PLJTF ' IGCB, ID, -.LOG, I03ur ,20) 

Cl  Co  CALL.  SETA?; IGCB, 2.0} 

0169  CALL  VIEW?( IGCB, 0., 135., 0., 100.) 

0170  CALL  WIKDW(  IGCB.O. ,  150. ,0. ,  1-00. ) 

0171  CALL  CSIZE(IGC3,3.) 

0172  CALL  FXD( IGC3,0) 

0173  CALL  PEN( IGCB, 1 ) 

0174  C 

0175  C  SCALE  DATA 
0176  C 

0177  IEXP=0 

0178  JEXP=0 

0179  IF ( I AXES . ME . 0 . AND . I AXES . NE . 3 ) GOTO  215 

0180  CALL  SCALE  ( X , N? ,  XMIM ,  XMAX ,  XTIC ,  IZEF.O .  I  EXP ,  LU ,  ISCAL ) 

0181  215  IF( I AXES. NE.O. AND. IAXES.NE. 2) GOTO  220 

0182  CALL  SCALE  ( Y .  NP ,  YMIN ,  YKAX ,  YTIC ,  IZERO ,  JEXP ,  LU ,  ISCAL ) 

0183  220  IF( IAXES.NE. 1 . AND. I AXES. NE. 2) GOTO  230 

0184  CALL  LSCAL( X , NP , XMIN , XMAX , XTIC , LU . ISCAL , LERR ) 

0185  IF ( LERR. EG. 1) GOTO  490 

01 80  E0225  K= 1 , NP 

0187  X(K)sALOCTCXCK)) 

0188  IF(DELX(K).SO.O.)GOTO  225 

0 1 89  DELX (K)  = ALOGT ( AB3 ( DELX (K) ) ) 

0190  225  CONTINUE 

0191  230  IF( IAXES.NE. 1 .AND. IAXES.NE. 3)GOTO  240 

0192  CALL  LSCAL(  Y, NP ,  YMIN , YKAX ,  YTIC ,  L'J ,  ISCAL ,  LERR ) 

0193  IF ( LERR. EQ.1) GOTO  490 

0194  D0235  Ks 1 , NP 

0195  Y(K)=ALCGT(Y(K) ) 

0196  IF(DELY(K) .EQ.O. )GOTO  235 

0197  DEL Y ( K ) =  ALOGT ( ABS( DELY 1 K ) ) ) 

0198  235  CONTINUE 

0199  C 

0200  C  WRITE  AXES  LABELS  AND  TITLE 
0201  C 

0202  240  CALL  MCVE( IGCB , 35 . , 1 . ) 

0203  IF( IEXP.NE.O)  GOTO  250 

0204  CALL  LAEEL( IGCB) 

0205  WRITE(Li;G,245)IAX 

0206  245  FORMAT(15A2) 

0207  GOTO  260 

0208  250  CALL  LABEL ( IGCB) 

0209  'WRITE;  LUG,  255)  IAX.IEXP 

0210  255  F0RNATO5A2,"  X  10**",  12) 

0211  260  CALL  KOVE( IGCB , 3 . . 1 2. ) 

0212  CALL  LDIR( IGCB, +1.57) 

0213  IF(-JEXP.NE.O)  GOTO  265 

0214  CALL  LABEL( IGCB) 

0215  WRITE( LUG, 245 : I AY 

0216  GOTO  270 

0217  265  CALL  LABEL; IGCB) 

0218  WRITE; LUG, 255) IAY, JEXP 

0219  270  CALL  MOVE  ( IGCB .  40 . .  90 . ) 

022C  CALL  LDIR'IGCB.C.) 
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0CC5 

LPLOT  GjVO  AM  MON. ,  26  JAN.,  1?Sl 

C221 

CALL  LABEL: IGCB } 

222- 

WRITE: LUG , 276 ) IKED 

0223 

1  i  y 

FORMAT: 20A2) 

0224 

( ' 

w 

022b 

■_  Dr. 

AW  AXES 

02-'; 

c 

0221 

IF ( IAXES . ME. 1 . AND. I AXES . NE . 2 ) GOTO  280 

■j  226 

CALL  LCOAXC 1 .XKIN,XMAX,XTIC,LUG) 

0229 

2*30 

IF ( I AXES . NE . 1 . AND . IAXES . NE . 3 ) GOTO  285 

0230 

CALL  LOGAX ( 2 , YMIN , YMAX , YTIC , LUG ) 

0231 

2g- 

IF ( IAXES . NE. 0 . AND . IAXES . ME . 3 ) GOTO  2<>0 

0232 

CALL  LINAX( 1 ,XMIN,XXAX,XTIC,LUG) 

023? 

290 

IF ( I AXES . NE .0. AND. I AXES . NE . 2 ) GOTO  2^5 

0234 

CALL  linax(2,ymin,ymax,ytic,lug) 

0235 

295 

CALL  VI£WP( IGC3, 17. , 120. , 10. ,80. ) 

0236 

CALL  WINDW( IGCB , XMIN , XMAX , YMIN , YMAX 1 

0237 

r~ 

0238 

C  PLOT  DATA  POINTS 

0239 

c 

^24C 

3C0 

MM.=0 

0241 

IF<  JJ.EQ.2):v!=1 

0242 

CALL  LINE: IGCB, MM) 

0243 

X310  X=1,NP 

0244 

IF(K.EQ. 1 .OR. JJ.EC.3)  GOTO  305 

0245 

CALL  DRAW( IGCB, X(K),Y(K)) 

024o 

C 

MOVE  TO  (X, Y) .CENTER  CHAR., PLOT  CHAR. 

0247 

305 

CALL  MOVE: IGCB, X'K) ,Y(K) ) 

0243 

IF(JJ.EQ.I)  GOTO  ’10 

0249 

CALL  C?LOT( IGCB,-HW,-HH,-2) 

0250 

CALL  LABEL (IGCB) 

0251 

WRITE ( LUG , 3 1 5 ) JCHAR( K K 5 

0232 

CALL  M.OVE(  IGC3 ,  X '  K ) ,  Y(  K ) ) 

0253 

310 

CONTINUE 

025* 

315 

FCRMAT( 1A1) 

0255 

CALL  PEN UP ( IGCB) 

0256 

CALL  LINE( IGCB, 0) 

0257 

write: LU, 320) 

0258 

320 

FORMAT:"  LEAST  SQUARES  LINE  ?",/.22X.' 

0259 

1  22X, "YES  TYPE  1") 

0260 

read: lu,*'  iline 

0261 

WRITE(LU, 325) 

0262 

325 

FORMATC  lRRQR  BARS  ?",/, 

0263 

1  "  NONE  TYPE 

0264 

’  "  X  TYPE 

0265 

3  "  Y  TYPE 

02o6 

4  "  X  4  Y  TYPE 

0267 

read: lu,*) ibar 

0268 

if: ILINE. EO.O)  GOTO  360 

0269 

d 

0270 

Z  CALCULATE  START  1 

027  ■ 

0272 

CALL  LSREGfX, 

0273 

XST=XMIN 

027 J 

YSTsYINT 

0^  j  } 

y3T=:glc?e*>: 
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027o  ;?0  IF( YINT.LZ. WAY.'  GO  PC  335 

02”7  XST= ( YMAX-YIN T ) /3LCPE 

0278  Y3T=YMAX 

0279  335  IF(YINT.GE.YMIN)  GO  TO  340 

0280  XST=(YMIN-YINT)  /SLOPE 

0281  YST=YKIN 

0282  340  XEND=XKAX 

0283  YEND= ( 5L0PE*XMAX ) + YINT 

0234  IF( YEND.LE.YMAX)  GO  TO  345 

0285  XEND= ( YMAX- YINT ) /SLOPE 

0256  YENDsYMAX 

0267  345  IF(YEND.GE.YMIN)  GO  TO  350 

0268  YENDsYKIN 

0289  XEND= ( YKJN- YINT ) /SLOPE 

02  9Q  C 

029 1i  C  PLOT  LEAST  SQUARES  Lli.E 
029 2  C 

0293’  350  IFULINE.EC.01  GOTO  760 

0294  CALL  MOVE ( IGCB , XST , Y3T ) 

02°5  CALL  DRAW ( IGCE ,  XEND ,  YE.ND ) 

0290  CALL  PENUP(IGCB) 

0297  IF( IAXES.NE.O)GOTO  360  1 

0293  LUX=LU 

0299  IF  CD.  EC.  2)  LUX=o 

0300  SLOPE=SLO?E» ( 1 0 . ** ( ( JEXP-IEXP ) * ( - 1 ) ) ) 

030 1  YINT= YINT* ( 1 0 . •* ( ( JEXP ) * ( - 1 ) ) ) 

0 202  S1aS1*f 10.**((JEXP-IEXP)*(-1))) 

0303  S2=S2*( 10.**( ( JEXP)*( -1 ) ) ) 

0304  WRITE ( LUX , 355 ) SLOPE , S 1 , YINT , £2 

0305  355  FORMAT ( "  SLOPE="G12.5",  *-"012.6/"  Y-ir.t="G12.5",  *-”012.6) 
0306  C 

0307  C  COMPUTE  AND  PLOT  ERROR  BARS 

0308  C 

0309  360  IF( I3AR.EQ.0)  GOTO  465 

0310  C  SCALE  DEVIATIONS  BY  VALUE  OF  EXPONENT 

0311  DO390  Ksl.NP 

0312  IF ( IEXP ) 365 , 375 , 370 

0313  365  CELX(K)=ABS(DELXf  )<)/(  10.*#IABS(IEX?) )) 

0314  GOTO  375 

0315  370  DELX(K)=ABS(DELX(K)*( 10.**IABS( IEXP) ) / 

0316  375  IF (JEXP) 380, 390, 335 

0517  380  DELY(K)=AE5(DELY(K)/( 10.*#IAES( JEXP) ) ) 

0313  GOTO  3Q0 

0319  385  DELY(K) =ABS(DELY(K)*( 10.**IABS( JEXP ) ) ) 

0320  390  CONTINUE 

0321  NOBARsO 

0322  0  DON'T  DRAW  EPROR  BARS  LESS  THAN  1.05?  X  OR  1.25%  Y  OF  LINEAR  AXIS 

0323  EARX=  ( XT  1AX-XKIN )  *0 . 0 1 05 

0324  3 ARY  = ( YMAX- YKIN ) *0 . 0 1 25 

0325  X455  K=1,NP 

0326  IF  ( IAXES.NE.  1  .AND.  I  AXES.  N'E.  2  5  GOTO  395 

0327  C  XN’T  DRAW  ERROR  BARS  LESS  THAN  13.5*  X  OR  18%  Y  OF  LOO  POINT 

0323  3ARX=(ABS( 10.**X(X) ) )»0. *35 

0329  XBAR1=AL0GT(ABS(  ( 1Q.#*X(K)  )-•'  13.«EELX(K; ) '  > 

0330  XBAR2=AL0CT( ABS( ( iQ . **X' K) ) ♦( 1C . *«D£LX ( K  V ) i 
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w  ^  —  1 

3 '.-6 

IF I  AXES .  ME . 1 .  AND .  I AXES .  NE .  3 )  GOTO  400 

c  ■?‘?2 

3ARY=( A3S( 10.**YCO))*0.l8 

5^23 

Yh An 1 = ALOGK ABS( ( 1 0 . *#Y( K) ) +( 1 0 . **DELY( K ) ) ) ) 

c:3* 

YEAR2= -LOGT(  ASS(  ( 1Q.**Y(K) )-(  !0.#»DELY(i<) ) ) ) 

0';6 

400 

IF ( I AXES . NE . 0 . AND . I AXES . NE . ? ) GOTO  405 

3326 

XEAR1=X(K)-DELX(K) 

0237 

XE AR2= X ( K) +DELX ( K ) 

0328 

*05 

IF { I AXES. NE. 0. AND. I AXES. NE. 2) GOTO  410 

3239 

YBAR1 =Y(K)+DELY(K) 

32-0 

YBAR2= Y ( K ) -DEL Y ( K ) 

32*1 

410 

3GTO( 420 , 440 , 420 )  IEAR 

02*2 

115 

GOTO  465 

0343 

C 

DRAW  X  ERROR  BARS 

334* 

-20 

IF( I AXES . NE . 1 .AND. I AXES. NE. 2) GOTO  425 

■3345 

IF( 1 0. *»DELX (K).LT.3ARX) NOBAR: 1 

034o 

IF( 10.#*DELX(K) .LT.3ARX1GOTO  435 

0347 

GOTO  430 

0348 

*25 

IF(DSLX(K).LT.BARX)NOBAR:1 

03*** 

IF(DELX(iC).LT.BARX)GOTO  435 

0350 

430 

GALL  MOVE( IGC3.X(K) ,Y(K) ) 

0381 

CALL  CPLOTCIGCB, -0.5,0. ,-2) 

0352 

CALL  DRAW( IGC3 , XBAR 1 , Y( K) ) 

03>3 

CALL  CPL0T(IGCB,0..-0.3,-2) 

0354 

CALL  CPL0T(IGCo,0.,0.6,-1) 

0335 

CALL  MOVE( IGC3,X(K) ,Y(K) ) 

O33o 

CALL  CPLOT( IGCB, 0.5,0. ,-2) 

0357 

CALL  DRAW( IGC3 , XBAR2 , Y( X ) ) 

038<5 

CALL  CPLOT ( IGC3 , 0 . , -0 . 3  > -2 ) 

0369 

CALL  CPLOT(IGCE,0.,0.6,-1) 

0360 

435 

IF( IBAR.EO. 1 )  GOTO  455 

03oi 

C 

DRAW  Y  ERROR  3ARS 

0362 

440 

IF( IAXES.NE. 1 .AND. IAXES.NE. 2 /GOTO  445 

0363 

IF ( 1 0 . »*DEL Y  C  K ) . LT . BARY ) NOBAR: 1 

0264 

IF(10.**DELY(K).LT.3ARY)GCTO  455 

03c5 

GOTO  450 

0366 

445 

IFOELY(K) .LT.3ARY)  NOBAR: 1 

0367 

IF(DELY(X) .LT.BARY)  GOTO  455 

0368 

450 

CALL  MOVEf ICCB,X(X) , Y(K) ) 

O'.t'J 

CALL  CPLOT ( IGCB . 0 . , 0 . 3 , -2 ) 

0370 

CALL  DRAW( IGCB.X(K) , YBAR1 ) 

0371 

CALL  CPLOT( IGCB, -0.3,0., -2) 

0372 

CALL  CPLOT{ IGCB, 0.6,0. ,-1 ) 

0372 

CALL  MOVE( IGCE.X(K) ,Y(K) ) 

0374 

CALL  CPLOT ( IGCB, 0. ,-0.3, -2) 

0378 

CALL  DRAW( IGCB , X( K ) , YBAR2 ) 

C37g 

CALL  CPLOT ( IGCB , -0 . 3 , 0 . , - 2 ) 

0377 

CALL  CPL0T( IGCB, 0.6,0. ,-1 ) 

0370 

455 

CONTINUE 

0379 

CALL  PENUP(IGCB) 

03-50 

IF ' N03AR . EQ .  1 )  WRITE (LU,4oO) 

0381 

*o0 

FORMAT ("  ***  SOME  ERROR  BARS  WERE  TOO  SMALL  TO  PLOT  ***") 

0282 

at5 

WRITE!L'J,470) 

0383 

470 

FORMAT ("  TO  EXIT  TYPE  0",/, 

0384 

1  "  IF  YOU  WANT  ANOTHER  PLOT  ON  SAME  AXES"./, 

0:35 

2  ”  AND  ALL  X  i  Y  VALUES  ARE  WITHIN  THE",/, 
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0336 

3  "  SCALES  OF  THE  FIRST  PLOT 

0367 

READ ( LU , * ) IDONE 

0363 

IF(KK.EQ.6)  GOTO  490 

0369 

C 

SELECT  PEN  COLOR  IF  USING 

PLOTTER 

0390 

IF(ID.NE.2.0R.IDGNE.NE. 1 ) 

GOTO  435 

0391 

475 

WRITE(LU,430) 

0392 

480 

FORMAT ("  ***  PEN  COLOR  SELECTION  ***", 

0393 

1  "  BLACK 

TYPE  1",/, 

0394 

2  "  RED 

TYPE  2",/, 

0395 

3  "  GREEN 

TYPE  3",/, 

0396 

4  "  BLUE 

TYPE  4") 

0397 

READ(LU,«)  IPEN 

0398 

IF(IPEN.LT. 1 .0R.IPEN.GT.4) 

GOTO  475 

0399 

CALL  PEN(IGC3,IPEN) 

0400 

435 

IF( IDONE. EQ. 1 )  GOTO  20 

0401 

490 

CALL  PEN(IGCB.O) 

0402 

CALL  PLOTR(IGCB,ID,0) 

0403 

END 
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0001  F'li.4  ,L 

0002  SUBROUTINE  SCALE  ( A, HP,  AMIN.AKAX,  TIC,  IZERO,  IEXP.LU.  ISCAL i 

000-  C  THIS  ROUTINE  COMPUTES  MAX  &  MIN  VALUES .SCALES  DATA  TO 

0004  C  E  FORMAT  AND  DETERMINES  THE  NUMBER  OF  TICK  MARKS  PER  AXIS 


0005 

C 

A 

-  ARRAY  TO  3E  SCALED 

0006 

C 

NP 

-  NUHEER 

OF  POINTS  IN  ARRAY  A 

0007 

C 

AMIN 

-  MINIMUM 

VALUE  OF  A 

0008 

C 

AMAX 

-  MAXIMUM 

VALUE  OF  A 

0009 

c 

TIC 

-  NUMBER 

OF  TICK  MARKS  ON  AXIS 

0010  C  IZZRO  -  SET  TO  0  TO  FORCE  ORIGIN  TO  ZERO, NORMALLY  =  1 

0011  C  I EXP  -  EXPONENT  OF  BASE  10  TO  WHICH  A  IS  RAISED 

0012  C  LU  -  LOGICAL  UNIT  NUMBER  OF  TERMINAL 

0013  C  ISCAL  -  SCALING,  C  FOR  AUTOMATIC,  1  FOR  MANUAL 

0014  DIMENSION  A(256) 

0015  ISML=0 

0016  C  FIND  MAX  AND  MIN  VALUES  OF  A 

0017  IF( ISCAL, EQ. 1 ) GOTO  25 

00  Id  AMAX=A( 1 ) 

0019  AMINsA(l) 

0020  DOIO  1=1, NP 

0021  IF(A( I) .GT.AMAX) AMAX=A( I) 

0022  IF( A( I) ,LT. AMIN)AMIN=A( I) 

0023  10  CONTINUE 

0024  IF (AMAX.NE. AMIN) GOTO  15 

0025  AMAX= AMAX-t-0 .  0 1 

0026  AMIN= AMIN-0. 01 

0027  15  IF ( IZERO. NE.O. OR. AMIN. GE.O)  GOTO  25 

0028  WRITE (LU, 20) 

0029  20  FORMAT ( "  **  DATA  HAS  NEGATIVE  VALUES  ORIGIN  CAN'T  BE  ZERO 

0030  IZERO= 1 

0031  25  IF ( IZERO . EC . 0 ) AMIN=0 . 

0032  IEXP=0 

0033  DIF=ABS(  AMAX-AMIi.) 

0034  C  IF  DIF  IS  LESS  THAN  1/2  MAGNIFY  SCALE 

0035  IF(DIF.GT.0.5)GOTC  40 

00 3o  ISML=1 

0037  30  D035  Xs1,X P 

0033  35  A(K)=A(K)*10. 

0039  AMIN=AMIN*10. 

0040  AMAX=AMAX*10. 

0041  IEXPsIEXP-t-1 

0042  DIF=DIF*10. 

0043  IF(DIF.LT. 1 . )GOTO  30 

0044  GOTO  DO 

‘0045  C  SCALE  DATA  TO  E  FORMAT  AND  SAVE  EXPONENT  I  EX?) 

0040  40  :FCABS(AMAX).GE.1.0F.ABS(AMIN).GE.1)  GOTO  50 

0047  D045  K=1,NP 

0043  45  A(K)=A(K)*10. 

0049  AMIN=AKIN*10. 

0050  AMAX=AKAX*1C.  \ 

0051  IEXP=IEXP-t-1  \ 

0052  GOTO  40 

0053  50  IFr ASS( AKAX) .L£. 1Q0C. AND. A££( AMIN) .LE. 1000)  XTC  cO 
0C54  D055  K=1,NP 

0053  55  A(K;=A(K)/10. 
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0056 

AHIN=AhIN/10. 

0057 

AMAX:AMAX/10. 

0056 

I£XP=I£XP-1 

0059 

GOTO  50 

OOoO 

c 

DETERMINE  INTERVAL  FACTOR 

0061 

50 

DIF= ABS( AMAX-AMIN ) 

0062 

IF(DIF.GT.5.)GOTO  65 

0063 

FACT3=1 .0 

00o4 

GOTO  100 

0065 

65 

IF(DIF.GT. 10. )GOTO  70 

0066 

FACTR=2 . 0 

0067 

GOTO  100 

0063 

70 

IF(DIF.GT.25.)GOTO  75 

00o9 

FACTR=5.0 

0070 

GOTO  100 

0071 

75 

IF ( DIF. GT. 50) GOTO  80 

0072 

FACTR=10.0 

0073 

GOTO  100 

0074 

60 

IF(DIF.GT. 125)  GOTO  85 

0075 

FACTR=25 . 

0076 

GOTO  100 

0077 

85 

IF ( DIF. GT. 250) GOTO  90 

0078 

FACTR=50. 

0079 

GOTO  100 

0080 

90 

IF(DIF.GT.500)GOT0  95 

0081 

FACTR=100. 

0082 

GOTO  100 

0083 

95 

FACIR=200 . 

0064 

C  HANDLE  ROUND-OFF  ERROR 

0085 

100 

0NE=0. 99999 

0086 

IF( ISML.EQ. 1 )0NE=0.0 

0037 

IF( AMIN.GE.O. )GOTO  105 

0088 

c  •***  for  NEGATIVE  NUMBERS 

0089 

TRUNK= ( AMIN/FACTR ) -ONE 

0090 

AMIN=AINT( TRUNK) *FACTR 

0091 

IF(AMAX.GE.O)  GOTO  110 

0092 

TRUNK=AMAX/FACTR 

0093 

AMAX= AINT ( TRUNK ) *FACTR 

0094 

GOTO  115 

0095 

C  ****  FOR  POSITIVE  NUK3ERS 

0096 

105 

TRUNK= AMIN/FACTR 

0097 

AMIN= AIMT( TRUNK ) *FACTR 

0098 

110 

TRUNK: ( AMAX/FACTR ) +CNE 

0099 

AKAXs AINT ( TRUNK ) *FACTR 

0100 

C  a#***#####*######*###*##***# 

0101 

115 

DIF: ABS( AMAX-AMIN) 

0102 

C 

DETERMINE  NUMBER  OF  TICK 

0103 

IF(DIF.NE. 1 . )GOTO  120 

0104 

TIC: 1 . 

0105 

GOTO  150 

0106 

120 

IF(DIF.GT. 10)GOTO  125 

0107 

DIF:DIF»10. 

0108 

GOTO  120 

0109 

125 

IF( DIF.LE. 100)GOTC  130 

0110 

DIF: DIF/ 10. 

JAN..  1931 


«*** 


#«** 
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0111 

GOTO  125 

0112 

130 

nc=3. 

0113 

IF( DIF.EQ.  100.  )G(JTO 

145 

0114 

:f(dif.eq.8o.)goto 

135 

0115 

IF(DIF.GE.60. )GOTO 

150 

0 1 T6 

IF( DIF.EQ. 50. )GOTG 

145 

0117 

IF( DIF.EQ. 40. )GOTO 

135 

0118 

IF( DIF.EQ. 30. )GOTO 

150 

0119 

IF (DIF.EQ. 25. )GOTO 

145 

0120 

IF (DIF. EC. 20.) GOTO 

140 

0121 

IF( DIF.EQ. 15. )GOTO 

150 

0122 

IF(DIF. EQ. 12. 5)GCTC 

1  1  45 

0123 

135 

TIC=4. 

0124 

GOTO  150 

0125 

140 

TIC=2. 

0126 

GOTO  150 

0127 

145 

TIC=5. 

0128 

150 

CONTINUE 

0129 

RETURN 

0130 

END 

FTN4  CCPPILER:  HPQ2060-16092  REV.  2001  (791101) 
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0001 
0002 
GOO  j 
^004 

G005 

jCO'j 

0007 
OOOd 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0015 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
3043 
0044 
0045 
00 4 o 
0047 
004o 
0049 
0050 
005 1 
0052 
0053 
0054 
0055 


FTM4 , L 


r 


r* 


SUBROUTINE  LSCAL  ( A , MP , PMIN , ?MAX , TIC , LU , ISCAL , LERR ) 

LOGARITHMIC  SCALING  ROUTINE  TO  COMPUTE  MAX  &  MIN  TO  LOG  BASE  10 

AND  DETERMINE  THE  NUMBER  OF  TICK  MARKS  PER  AXIS 

A  -  aRRAY  TO  BE  SCALED 

N?  -  NUMBER  OF  POINTS  IN  ARRAY  X 

PMIN  -  POWER  TO  BASE  10  OF  MINIMUM  VALUE  OF  ARRAY  A 

PMAX  -  POWER  TO  BASE  10  OF  MAXIMUM  VALUE  OF  ARRAY  A 

TIC  -  NUMBER  OF  TICK  MARKS  ON  AXIS 

lU  -  LOGICAL  UNIT  NUMBER  FOR  ERROR  OUTPUT 

LSCAL  -  SCALING,  0  FOR  AUTOMATIC,  1  FOR  MANUAL 

uERri  -  LOG  ERROR,  CAN’T  TAKE  LOG  OF  A  NON-POSITIVE  NUMBER 

DIMENSION  A( 256 ) 


r 


C  ILST  ARRAY  FOR  MON-POSITIVE  VALUES 


LERR=0 
D01D  1=1, NP 
IF(Ad).GT.O)  GOTO  10 
LERR=1 

10  CONTINUE 

IF  ( LERR .  EQ .  1 )  WRITE(  LU ,  1 5 ) 

IF ( LERR. EC. 1) GOTO  90 
IF ( ISCAL . EC . 1 )GOTO  25 

15  FORMAT ("  STOP  YOU  CAN'T  USE  LOGARITHMIC  SCALING  ROUTINE  (LSCAL)" 

-/"  FOR  NON-POSITIVE  VALUES,  USE  A  LINEAR  SCALING  ROUTINE") 

C 

C  DETERMINE  VALUES  OF  MAX  AND  MIN 
C 

AMIN=A( 1 ) 

AMAXsAMIN 
D020  1=2, NP 

IF (AMIN . GT . A( I ) ) AMIN=A( I ) 

IF ( AMAX . LT. A( I ) ) AMAX= A( I ) 

20  CONTINUE 

GOTO  30 

25  AMINsPtCTN 

AMAXsPMAX 

C 

C  COMPUTE  LOG  OF  MIN  AND  RAX  VALUES 


30  ?OWL=  ALOGT ( AMIN ) 

MINPrPOWL 
POWHs ALOGT (AMAX) 

MAXP=POWH 

IF ( ABS( FLOAT ( MAXP ) -POWH ) . LT . 1 . 0E-05 ) GOTO  35 
IF(POWH.GE.O.O)MAXP=MAXP+1 
35  IF ( ABS( FLOAT (M1NP) -FOWL ).LE.1.0E-05) GOTO  40 

IF ( FOWL. GE.O) GOTO  40 
MINP=MINP-1 
40  ?MIN=MINP 

PMAX=MAXP 

C  MAKE  DIFFERENCE  IN  EXPONENTS  OF  MAX  i  MIN  A  MULTIPLE  OF  1,2, 3, 4  OR  5 
IDIF= ABS ( PMAX-PMIN ) 

IF(IDIF.LE.6)G0T0  55 
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0056  IF ( IDIF . EQ . 7 . OR . IDIF . EQ . 9 ) PMAX=PMAX+ 1 

0057  IF ( IDIF . LE. 1 0 )GOTO  55 

0058  LDIF=15 

0059  45  IF (LDIF.GE. IDIF) GOTO  50 

0060  LDIF=LDIF+5 

0061  GOTO  45 

0062  50  PMAX=PMIN+LDIF 

0063  C 

0064  C  DETERMINE  NUMBER  OF  TICK  MARKS  PER  AXIS 
0065  C 

0066  55  IDIF  =  ABS ( PMAX-PMIN ) 

0067  IF(IDIF.NE. 1 )GOTO  60 

0068  TIC=1 

0069  GOTO  90 

0070  oO  IF(IDIF.GT. 10)GOT0  65 

0071  IDIF=IDIF*10 

0072  GOTO  60 

0073  65  IFdDIF.LE.  100) GOTO  70 

0074  IDIFsIDIF/10 

0075  GOTO  65 

0076  70  TIC=3 

0077  IF (IDIF. EQ. 100) GOTO  85 

0078  IF (IDIF. EQ. 80) GOTO  75 

0079  IF ( IDIF. GE. 60) GOTO  90 

0080  IF (IDIF. EQ. 50) GOTO  85 

0081  IF(IDIF.£Q.40)GOTO  75 

0082  IF(IDIF.EQ.30)G0T0  90 

0083  IF ( IDIF. EQ. 25) GOTO  85 

0084  IF (IDIF. EQ. 20) GOTO  80 

0085  IF( IDIF.EQ. 15)GOTO  90 

0086  75  TIC=4 

0087  GOTO  90 

0088  80  TIC=2 

0089  GOTO  90 

0090  85  TIC=5 

0091  90  CONTINUE 

0092  RETURN 

0093  END 

FTN4  COMPILER:  HP92060- 1 6092  REV.  2001  (791101) 

**  NO  WARNINGS  **  NO  ERRORS  «  PROGRAM  =  00497  COMMON  =  00000 
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PAGE  0001  FTN.  2:3o  PM  TUS.,  30  DEC.,  1980 
0001  FTN4.L 

0002  SUaROUTINE  LINAX( 1 AXIS, AMIN, AMAX, TIC, LUG) 

0003  C  LINEAR  AXIS  DRAWING  ROUTINE 
0004  C 

0005  C  I AXIS  -  IrX  AXIS  ,  2=Y  AXIS 

0006  C  AMIN  -  MINIMUM  VALUE  OF  AXIS 

0007  C  AMAX  -  MAXIMUM  VALUE  OF  AXIS 

0008  C  TIC  -  NUMBER  OF  TICK  MARKS  ALONG  AXIS 

0009  C  LUG  -  LOGICAL  UNIT  NUMBER  FOR  GRAPHICS  OUTPUT 

0010  C 

0011  COMMON  IGCB092) 

0012  C  SET  AXIS  LENGTH  IN  WORLD  COORDINATE  SYSTEM  (WCS) 
0013  IF( IAXIS.SQ. 1 )ALEN= 114. 

0014  IF ( IAXIS . £Q . 2 ) ALEN=70 . 

0015  C  DEFINE  ORIGIN 
0016  X0= 1 9 • 

0017  Y0=10. 

0018  TICM=AMIN 

0019  CALL  MOVE ( IGCB , XO , YO ) 

0020  CALL  CSIZE(IGCB,3.) 

0021  IF( IAXIS.EQ.2)G0T0  100 

0022  C 

0023  C  DhAW  X  AXIS 
0024  C 

0025  CALL  DRAW( IGCB , ALEN+XO , YO ) 

0026  C  LABEL  X  ORIGIN 

0027  CALL  MOVE ( IGCB ,XO,YO) 

0028  CALL  MOVEI( IGCB,-12. ,-4.5) 

0029  CALL  LABEL ( IGC3) 

0030  WRITE (LUG, 30 )TICM 

0031  30  F0RMAT(F8.Q) 

0032  C  DRAW  X  TICK  MARKS 
0033  DO  50  K=1,TIC 

0034  T ICK=  ALEN* ( FLOAT ( K ) /T IC ) 

0035  CALL  MOVE ( IGCB, TICK+XO,YO) 

0036  CALL  DRAW( IGCB,TICK+XO, YQ+2.0) 

0037  C  LABEL  X  TICK  MARKS 

0038  TICM=TICM*( ( AMAX-AMIN) /TIC) 

0039  CALL  MOVE( IGCB,TICK+X0, YO) 

0040  CALL  MOVEI ( IGCB , - 1 2 . , -4 . 5 ) 

0041  CALL  LABEL (IGCB) 

0042  WRITE (LUG, 30 )TICM 

0043  50  CONTINUE 

0044  GOTO  200 

0045  C 

0046  C  DRAW  Y  AXIS 
0047  C 

0048  130  CALL  DRAWUaCB,XO,ALE!'H.YO) 

0049  C  LABEL  Y  ORIGIN 

0050  CALL  MOVE ( IGCB, XO,YO) 

0051  CALL  MQVEK IGCB, -17. ,-0.8) 

0052  CALL  LABEL(IGCB) 

0053  WRITE(LUG,30)TICM 

0054  C  DRAW  Y  TICK  MARKS 
0055  DO  150  Ksl.TIC 
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PAGE  0002  LIN AX  2:36PM  TUE. ,  30  DEC.,  I960 
0056  TICK= ALEN* ( FLOAT ( K ) /TIC ) 

0057  C  LABEL  Y  TICK  MARKS 

0056  TICMsTICM+((AMAX-AMIN)/TIC) 

0059  CALL  MOVE ( IGCB , XO , TICK+YO ) 

0060  CALL  DRAW( IGCB , XO+2 . 5 , TICK+YO ) 

0061  CALL  MQVE( IGCB, XO, TICK+YO) 

0062  CALL  MOVEI(IGCB,-17.,-0.8) 

0063  CALL  LABEL (IGCB) 

0064  WRITE(LUG,30)TICM 

0065  150  CONTINUE 

0066  200  CALL  PENUP(IGCB) 

0067  RETURN 

0068  END 

FTN4  COMPILER :  HP92060-16092  REV.  2001  (791101) 

»*  NO  WARNINGS  **  NO  ERRORS  **  PROGRAM  s  00342  COfflON  s  00192 
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PAGE  0001 

FTN.  2:23  PM  TUE. ,  2  DEC.,  1980 

0001 

FTN4, 

* 

0002 

SUBROUTINE  LOGAX( I AXIS, PMIN, PMAX, TIC, LUG) 

0003 

C  LOGRATHMETIC  AXIS  DRAWING  ROUTINE 

0004 

C 

0005 

C  IAXIS  -  1=X  AXIS  ,  2=Y  AXIS 

0006 

C  PMIN  -  POWER  TO  BASE  10  OF  MINIMUM  VALUE  OF  AXIS 

0007 

C  PMAX  -  POWER  TO  BASE  10  OF  MAXIMUM  VALUE  OF  AXIS 

; 

0008 

C  TIC 

-  NUMBER  OF  TICK  MARKS  ALONG  AXIS 

i 

0009 

C  LUG 

-  LOGICAL  UNIT  NUMBER  FOR  GRAPHICS  OUTPUT 

0010 

C 

i 

0011 

COWON  IGCB(192) 

0012 

C  SET  AXIS  LENGTH  IN  WORLD  COORDINATE  SYSTEM  (WCS) 

0013 

IF(IAXIS.EQ. 1 )ALEN=114. 

0014 

IF( IAXIS . EQ . 2 ) ALEN=70 . 

0015 

C  DEFINE  ORIGIN 

i 

0016 

X0=19. 

i 

0017 

Y0=10. 

0018 

IEXP=PMIN 

0019 

CALL  CSIZE(IGCB,3.) 

0020 

CALL  MOVE( IGCB , XO , YO ) 

0021 

IF(IAXIS.EQ.2)G0T0  25 

0022 

C 

0023 

C  DRAW  X  AXIS 

0024 

C 

0025 

CALL  DRAW ( IGCB , ALEN+XO , YO ) 

0026 

C  LABEL  X  ORIGIN 

i 

0027 

CALL  MOVE ( IGCB, XO,YO) 

0028 

CALL  MOVEI ( IGCB , -3 . 7 , -4 . 5 ) 

0029 

CALL  LABEL (IGCB) 

0030 

WRITE(LUG.IO) 

0031 

10 

FORMATC'IO") 

0032 

CALL  CSIZE(IGCB,2.) 

0033 

CALL  MOVEI( IGCB, 2. 5, 2. ) 

0034 

CALL  LABEL (IGCB) 

0035 

WRITE(LUG, 15)IEXP 

0036 

15 

F0RMAT(I3) 

0037 

CALL  CSIZE( IGCB , 3 • ) 

0038 

C  DRAW  X  TICK  MARKS 

0039 

D020  K= 1 , TIC 

0040 

TICK=ALEN*(FLOAT( K) /TIC ) 

0041 

CALL  MOVE( IGCB , TICK+XO , YO ) 

0042 

CALL  DRAW ( IGCB , T ICK+XO , YO+2 . 0 ) 

0043 

C  LABEL  X  TICK  MARKS 

0044 

CALL  MOVE ( IGCB, TICX+XO.YO) 

0045 

CALL  MOVEI (IGCB, -3. 7, -4. 5) 

l 

0046 

CALL  LABEL (IGCB) 

0047 

WRITE(LUG, 10) 

0048 

CALL  CSIZE(IGCB,2.) 

0049 

CALL  MOVEI( IGCB, 2. 5, 2. ) 

0050 

IEXP= IEXPf ( PMAX-PMIN ) /TIC 

0051 

CALL  LABEL (IGCB) 

0052 

WRITE( LUG ,15) IEXP 

0053 

CALL  CSIZE(IGCB,3.) 

0054 

20 

CONTINUE 

l 

0055 

GOTO  35 

l _ 
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PAGE  0002  LOGAX  2:23PM  TUE.f  2  DEC.,  1980 
0056  C 

0057  C  DRAW  Y  AXIS 
0058  C 

0059  25  CALL  DRAW(IGCB,XO,ALEN+YO) 

0060  C  LABEL  Y  ORIGIN 

0061  CALL  MOVE(IGCB,XO, YO) 

0062  CALL  MOVEI(IGCB,-9.5,-1.5) 

0063  CALL  LABEL (IGCB) 

0064  WRITE(LUG, 10) 

0065  CALL  CSIZE(IGCB,2.) 

0066  CALL  M0VEI(IGCB,2.5,2.) 

0067  CALL  LABEL(IGCB) 

0068  WRITE(LUG,15)PMIN 

0069  CALL  CSIZE(IGCB,3.) 

0070  C  DRAW  Y  TICK  MARKS 

0071  D030  K= 1 , TIC 

0072  TICK=ALEN*(FLOAT(K)/TIC) 

0073  C  LABEL  Y  TICK  MARKS 

0074  CALL  MOVE ( IGCB , XO , TICK+YO ) 

0075  CALL  DRAW( IGCB , XO+2 . 5 , TICK+YO ) 

0076  CALL  MOVE( IGCB , XO .TICK+YO ) 

0077  CALL  MOVEI ( IGCB , —9 .5,-1 .5) 

0078  CALL  LABEL (IGCB) 

0079  WRITE(LUG.IO) 

0080  CALL  CSIZE(IGCB,2.) 

0081  CALL  M0VEI(IGCB,2.5,2.) 

0082  IEXPs IEXP+( PMAX-PMIN ) /TIC 

0083  CALL  LABEL (IGCB) 

0084  WRITE(LUG,15)IEXP 

0085  CALL  CSIZE(IGCB,3.) 

0086  30  CONTINUE 

0087  35  CALL  PENUP(IGCB) 

0088  RETURN 

0089  END 

FTN4  COMPILER:  HP92060- 1 6092  REV.  2001  (791101) 

»*  NO  WARNINGS  **  NO  ERRORS  *«  PROGRAM  =  00439  CONMON  =  00192 
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2  DEC.,  1980 


PAGE  0001  FTN.  1:44PM  TUS. , 

0001  FTN4,L 

0002  SUBROUTINE  L3REG  ( X , Y , IFRM, ITO , SLOPE , YINT , SDSLGP , SDYINT ) 

0003  C  PROGRAM  CALCULATES  LEAST  SQUARES  REGRESSION 
0004  C  X  -  X  ARRAY 

0005  C  Y  -  Y  ARRAY  i 

0006  C  IFRM  -  STARTING  POINT  OF  INTERVAL 

0007  C  ITO  -  ENDING  POINT  OF  INTERVAL 

0008  C  SLOPE  -  RISE  OVER  RUN  OF  L.S.  LINE 

0009  C  YINT  -  Y  INTERCEPT  OF  L.S.  LINE 

0010  C  SDSLOP  -  STANDARD  DEVIATION  OF  SLOPE 

0011  C  SDYINT  -  STANDARD  DEVIATION  OF  Y  INTERCEPT 

0012  DIMENSION  X(25o) , Y( 256) 

0013  FN= ITO- IFRM*  1 

0014  TX=0 

0015  ZY=C 

0016  XY=0 

0017  WY=0 

0018  SY=0 

0019  DO  100  KsIFRM, ITO 

0020  SY=SY+(Y(K)**2) 

0021  WY=WY+(X(K)*Y(K) ) 

0022  XY=XY+X(K) 

0023  ZY=ZY+Y(K) 

0024  100  TX=TX+(X(K)**2) 

0025  TY=(ZY«2)/FN  !  \ 

0026  XZY=(XY*ZY)/FN 

0027  XYN=(WY-m)«*2 

0028  XYD=TX-( (XY**2)/FN) 

0029  SX=(XY**2)/FN 

0030  SDSLOP =SQRT( ( SY-TY-( XYN/XYD ))/(( FN-2 . ) * ( TX-SX ) ) ) 

0031  SDYINT=SQRT( (SDSL0P**2)*(TX/FN) ) 

0032  SLOPE= ( (FN*WY)-(XY*ZY) )/( (FN*TX)-(XY«*2) ) 

0033  YINTs(ZY-(SLOPE*XY))/FN 

0034  RETURN 

0035  END 

FTN4  COMPILER:  HP92060-16092  REV.  2001  (791101) 

**  NO  WARNINGS  »*  NO  ERRORS  *•  PROGRAM  s  00315  COMMON  s  00000 
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